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ABSTRACT 


This thesis describes the single link flexible beam testbed at the CLaMS laboratory 
in terms of its hardware, software and linear model and presents two controllers, each 
including a hub angle PD feedback compensator and one augmented by a second 
static gain full state feedback loop, based upon a synthesized strictly positive real 
(SPR) output, that increases specific flexible mode pole damping ratios w.r.t the 
PD only case and hence reduces unwanted residual oscillation effects. Restricting 
full state feedback gains so as to produce a SPR open loop transfer function ensures 
that the associated compensator has an infinite gain margin and a phase margin of 
at least [-90, 90] degrees. Both experimental and simulation data are evaluated 
in order to compare some different observer performances when applied to the real 
testbed and to the linear model when uncompensated flexible modes are included. 
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INTRODUCTION 


This thesis is organized into five chapters describing the hardware, software, mod- 
eling, control design and test results for the single link flexible manipulator system 
at the CLaMS laboratory. The testbed is designed so as to enable the study of the 
identification and control of a real flexible system and the concept is similar to that 
of the earlier work done by Cannon and Schmitz (see e.g., [3]). 

Chapter one details the hardware which is centered on an aluminum beam 
(a plate), driven by a motor at one end while the other end is free, which was 
designed to exhibit low frequency poorly damped flexible modes. The real time 
controlling microprocessor is a^O Mhz Inmos Transputer mounted on a motherboard 
attached to a VME backplane : analog to digital and digital to analog converter 
'cards mounted on the same backplane provide external signal interfaces. Currently 
a hub angle sensor and beam strain sensors are used while two tip displacement 
sensors are being installed, namely a light sensor and an induction sensor. 

The second chapter describes the main software functions which implement 
linear and non-linear controllers and run on two transputers : three controllers 
are designed in this thesis and each is characterized by a fixed interval discrete time 
linear observer and compensator pair. Also, beam parameter identification is carried 
out with the same software. Reference trajectories and controller parameters axe 
loaded from Matlab data files before each experiment and afterwards the resulting 
data is stored in a similar fashion for later analysis. Before each experiment the 
D/A and A/D converter DC offsets may be loaded from another Matlab data file 
or updated via a short (about 25 seconds) test routine. 

Chapter three outlines the linear modeling of the beam, from the governing 
pair of linearized partial differential equations to a finite normal mode approximation 
in state space form and the computation of the output sensor matrix for the hub 
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angle potentiometer and beam strain gages. This linear model is characterized by a 
collection of second order transfer functions (the beam modes) driven by a common 
input torque and each sensor can be viewed as a linear combination of the associated 
states. The rigid body mode contains a stable real pole and a pole at the origin whilst 
each of the others have a very poorly damped but stable complex pole pair with 
a unique oscillation frequency (these are called the flexible modes). The complete 
transfer function is thus the linear combination of the individual modal transfer 
functions each of which dominates close to its center frequency. A brief description 
of the identification process is also given : direct modal parameter identification is 
used. 

The fourth chapter describes two full state compensators designed to improve 
the beam tip displacement response and each requires an observer to provide the 
state estimate. The first controller designed is a hub angle PD compensator which 
is found to give very good results in terms of tip displacement step response time, 
overshoot and residual oscillations. The design method used is simple : initially ig- 
nore the flexible modes and select those PD compensators which move the two rigid 
body mode poles such that the second order closed loop system is, say, critically 
damped. Then restore the flexible mode dynamics and select a single PD compen- 
sator from the previous set that also increases the flexible mode pole damping ratios 
with particular emphasis placed on the lowest frequency mode. 

PD compensation is a form of static gain full state feedback that has the 
advantage of an infinite gain margin in the ideal case (i.e. perfect modeling) and 
this will remain true for the actual plant closed so long as the torque to estimated hub 
angle velocity transfer function remains strictly positive real (SPR) as is guaranteed 
by the Passivity Theorem. However observer, sensor and other errors will make the 
actual system more fragile. 

The second controller designed is constructed by adding another static gain 
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feedback loop designed so as to further increase the damping ratios of the flexible 
modes. For this loop the time domain definition of strict positive realness is em- 
ployed to synthesize a SPR output from the state estimate of the PD compensated 
plant that emphasizes selected flexible modes. By restricting the class of feedback 
gains in this fashion we can guarantee that static negative feedback applied to this 
SPR SISO system will have an infinite gain margin and a phase margin of at least 
[—90, 90] (again in the ideal case). 

The actual gains are selected through trial and error however there is a strong 
correlation between the improvement in individual flexible mode pole damping ratios 
and corresponding terms in the synthesized SPR output matrix parameterization. 
Two such SPR outputs are selected, one with a modal transformation applied to 
the PD compensated system so as to explicitly decouple the flexible modes, and the 
other with the original coupled PD compensated system. 

Several closed loop observer designs are selected, namely a LQE full state ob- 
server, a pole placement full state observer and a reduced order Luenberger observer. 

In the fifth and final chapter some simulations and experimental results are 
examined in order to evaluate the performance of the controller/observer combina- 
tions in step response tests. The distorting effects of the observers in conjunction 
with unmodeled beam modes on the synthesized SPR transfer function is exam- 
ined through simulations and the observer which most faithfully reconstructs this 
transfer function (and performs well when applied to the true plant) is selected. Ex- 
perimental bode plots are also examined since they give a good indication of flexible 
mode pole damping ratios. 


CHAPTER 1 
Hardware 


1.1 Overview 

Figure 1.1 is an outline of the physical single flexible link testbed the central 
component being the flexible beam itself which is chosen to provide readily measur- 
able flexibility in the horizontal plane that is independent of gravity (ideally). This 
is driven by a DC motor at one end about a vertical axis in order to excite the flexible 
modes, and the other end is free. The motor is energized by a power amplifier which 
itself is controlled by a single voltage signal; the power amplifier-motor combination 
provides an approximately linear voltage to beam hub torque relationship. The hub 
angle is converted to a voltage signal by a potentiometer while the beam strains 
at specific locations axe likewise converted by strain gages. The realtime computer 
is a networked set of four Inmos Transputer micro-processors which are connected 
to a PC-AT host for software development and user interface. These Transputers 
also connect to a VME signal bus on which resides the digital to analog (D/A) and 
analog to digital (A/D) circuit cards providing the motor drive (computer output) 
and sensor (computer input) signal interfaces. 

1.2 Selecting the Flexible Beam 

A problem with a single flexible plate (i.e., a thin beam) is that the unsup- 
ported end can sag due to the effect of gravity when the beam is flexed and this is a 
source of nonlinear behavior and an unwanted complication. To minimize this effect 
a high Youngs modulus to density ratio is desired in order to maximize strength 
and minimize weight and mainly for this reason aluminum was selected as the beam 
material over steel, brass or plastic (see table 1.1); this ratio is similar for both steel 


1 



2 



Figure 1.1: Testbed Physical System 



and aluminum however the latter is much lighter and hence higher acceleration may 
be achieved for a given applied hub torque. Furthermore we can expect sagging to 
increase with a decreasing beam width or height to length ratio. 

a '■* Another factor to consider is the oscillation frequencies of the beam which 
' depend upon the material properties and the geometry of the system. For both the 
fixed-free (i.e. h = hub inertia -»■ oo) and the pinned-free (i.e. h = 0) cases the 
oscillation frequencies for the linear model are of the form (see equation 3.8), 

/• - IV? <“> 

where, 

f n = n t i j modal frequency, n = 0,1,.. 

E = Youngs Modulus, 
p = Density, 

L = Beam Length, 
w = Becim Width. 
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Table 1.1: Candidate Beam Materials 


Material 

units 

Density 

Kg/m 3 

Youngs Mod. 
Kg/m 2 

wgm 

Aluminum, Alloy 6061. 

2.713e3 

6.895el0 

5.041e3 

Steel, Type 302. 

8.028e3 

1.957ell 

4.937e3 

Brass, CA 230. 

8.748e3 

1.172ell 

3.660e3 

Plexiglas, typical values. 

1.18 e3 

2.969e09 

1.586e3 


When the hub inertia lies somewhere between the two extremes given above 
the equation for /„ becomes more complex however the values for /„ will fall between 
those of the fixed-free and pinned-free cases. 

In order to minimize the computer speed requirement low oscillation frequen- 
cies are desired and hence we wish to minimize the width to length ratio and min- 
imize the Youngs modulus to density ratio, however this will have the effect of 
increasing sagging at the tip and hence a compromise is required. 

A beam length of approximately one meter was chosen to provide a relatively 
large structure that will not be significantly affected by multiple sensors (e.g. strain 
gages, hub angle sensor, position light sensors) and to minimize the effects of hub 
shaft friction. Several beams of aluminum, steel and plastic were examined with a 
height of ten centimeters approximately and a width such that the fourth frequency 
was close to 50 Hz and the aluminum beam was selected as it sagged the least and 
had a low weight. 

1.3 Actuator 

A DC motor and switch mode power amplifier were selected to drive the beam 
as this combination provides a good degree of inherent open loop linearity at a 
reasonable cost and is a common method of actuating robot joints. The motor is a 
PMI Disc DC motor which contains no rotor iron other than the central shaft (the 
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coils are compressed into a flat disc) and this provides four main benefits, 

1. no inter-coil iron implies no cogging torque, a common unwanted non-linearity 
in other DC motor types, 

2. the rotor mass is concentrated near the shaft axis and hence has a low inertia 
which yields high acceleration; in the complete setup the hub mounting is the 
main hub inertial component, 

3. the small rotor coils have a low inductance value which increases the ampli- 
fier/motor bandwidth, 

4. low coil inductance eliminates commutation brush sparking which reduces mo- 
tor generated electrical noise. 

The main disadvantage of this particular motor is the shaft friction torque due 
to brush/coil contact and the use of solid bearings. 

A switch mode power amplifier, as its name implies, operates as a switch 
to alternatively apply full motor rotor coil voltage in one direction and then the 
other at a high frequency (20 Khz nominal). This is obviously a highly nonlinear 
process however the switching process is controlled so that the motor current is 
nearly linearly proportional to the input voltage signal over smaller bandwidth (e.g. 
I Khz) which is more than sufficient for our purposes. The only reasons for the 
switching operation are amplifier design considerations eg. power efficiency and 
component cost. 

The motor/amplifier combination provides a low pass filter voltage to torque 
transfer function and the dominant pole location is determined by the rotor circuit 
resis tanc e and inductance, 

t = l/LR (1*2) 

r = the dominant pole time constant, 
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L = rotor circuit inductance, 

R = rotor circuit resistance. 

and is estimated to be 2 msec approximately which means a 500 Hz bandwidth. 

While the low rotor inductance is desirable in order to maximize this band- 
width it has the unfortunate effect of increasing the motor rms current which de- 
termines the heat dissipated. A small inductor was added into the rotor circuit to 
contain this effect doubling the overall inductance (this is included in the 500Hz 
bandwidth estimation). 

1.4 Sensors 
1.4.1 Hub Angle 

A potentiometer was chosen for sensing the hub angle because it is cheap 
and because the interface circuitry is very simple (assuming an A/D converter is 
availible). However other options exist which can provide better resolution and 
accuracy and less noise disturbance e.g., an optical angular encoder or a resolver 
and one of these will eventually replace the potentiometer. 

Being an analog sensor the potentiometer has infinite resolution but its accu- 
racy is limited by circuit and contact noise and any non-linearity in the resistive 
material; the former affects relative accuracy while the latter reduces absolute ac- 
curacy. Furthermore, since the A/D converters have a 12 bit word length and 
the potentiometer translates +/ — 180 degrees into to the A/D full scale voltage 
of +/ - 10 V the actual hub angle resolution is limited to 0.0879 degrees approxi- 
mately. When the beam is at rest this corresponds to a tip displacement resolution 
of approximately 1.68mm . 
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1.4.2 Beam Strain 

Thin film strain gages are used to sense beam longitudinal strain at four points 
along the beam and the particular strain gages selected have the following features, 

1. a small size to provide a local measurement, 

2. a temperature drift coefficient to offset aluminum heat expansion, 

3. capable of measuring large repeated strains. 

In each case the strain gages are mounted in pairs on either side of the beam 
so that the resistance of one increases when the beam bends while the resistance 
of the other decreases; each pair forms part of a wheatstone bridge circuit and this 
arrangement provides good immunity to unwanted side effects such as temperature 
dependancies or contact potentials: the rest of the bridge circuit is comprised of two 
high quality matched resistors. 

1.5 Computer 

1.5.1 Real Time Computer 

The real time computer is a network of four Inmos Transputers (see figure 

/ - i 

2.1, only two are in use at present) which may be readily expanded in the future 
as proccessing requirements dictate. The Transputer is a single chip microprocessor 
which is availible in several formats and we use the T800-20 model which runs 
at 20 Mhz and is rated at 20 MIPS and 2.8 Mflops peak instruction rates. Each 
Transputer contains the following, 

1. 32 bit CPU 

2. 64 bit FPU 

3. 4 Kbytes of feist static RAM 

4. four 20 Mbit/sec bidirectional serial links 
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5. 32 bit data and memory buses and a DRAM controller 


The Transputers are specifically designed for external parallel operation using 
the high speed serial links for interconnections and these links also connect the 
transputers to the host development computer and the signal interface data bus. 

The transputers are also designed to facilitate multi-tasking (each separate 
program is called a process and high speed switching is achieved between separate 
processes on the same Transputer) and this constitutes a second form of parallelism. 
Separate processes communicate with each other through channels, not shared mem- 
ory, and if the processes lie on different Transputers then these channels are mapped 
to the physical hardware links before runtime. In this way the code may be re- 
distributed across different nets by simply defining a new mapping operation and 
therefore the processes need not be rewritten. Thus there are two distinct networks, 
the real hardware form and the desired software form. Note that while any number 
of software channels (memory permitting) may exist between processes running on 
the same Transputer, this is not the case for inter-Transputer communication as 
each Transputer has exactly four serial links. 

A third form of parallelism exists since the transputer allows the CPU, FPU 
and link interfaces to operate in parallel which further improves the computers 
performance if software is written to take advantage of this feature. 

The high level parallel language Occam was developed to make optimal use 
of the Transputers hardware features and thus experimental control and data in- 
put/output algorithms are written in this language. The code used was originally 
written in MicroSoft Quick C (almost ANSI C compatible) for an PC-AT386 based 
system and now an Inmos ANSI C compiler is used for non-speed critical sections 
while the rest has been converted to Occam. 
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1.5.2 Development Host Computer 

The user interface to the Transputers is provided by a PC-AT286 micro com- 
puter and when developing the software the Occam or ANSI C compilers and linkers 
are downloaded to the transputers along with the source code which then performs 
the compilation, linking and network configuring. At runtime the Transputers are 
reset and then the executable code is downloaded and booted. 

1.5.3 Computer to Signal Interface" 

1.5.3. 1 Interface Card Bus 

A VMEbus was selected to interface the Transputers to the A/D and D/A 
cards because of its large data bandwidth and the the wide variety of compatible 
circuit cards that exist and also because it is the bus used by an adjacent laboratory 
(the CIRSSE Dual Robot Arm Testbed). This bus supports 8, 16 and 32 bit data 
words at a maximum speed 40 MHz and the backplane is mounted along with the 
circuit power supply in a Eurocard rack. 

1.5. 3. 2 A/D Converters 

The A/D card selected is a XVME-566 made by XYCOM and features 32 
single ended, or 16 differential, input channels with 12 bit resolution and 12 bit 
accuracy. One channel may be sampled at a time (the conversion time is 10 fxsec) 
with a maximum sample throughput rate of lOOKhz and any number of channels 
up to the maximum may be sampled in any sequence and the results stored in on 
board RAM. 

In the current configuration single ended bipolar voltage input (+/ — 10 V) is 
used, the programmable gain is set to unity and the data ram is memory mapped 
to one of the Transputer microprocessors via the VMEbus so that each sample is 
accessed by a direct memory read command. 
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The inputs may be alternatively configured as 16 differential channels which 
should improve the signal noise common mode rejection ratio, however tests showed 
the measured noise to be greatly reduced when operated in single ended mode; I do 
not know why this is the case. 

1. 5.3.3 D/A Converters 

The D/A card is a VMIVME-4100 made by VME Microsystems Int. Corp. 
and it contains eight 12 bit +/ - 10 V output channels each of which has a settling 
time of Sfisec (max) to +/ — 1/2 LSB. Like the A/D card, the data RAM is memory 
mapped via the VMEbus to a Transputer and may be written to directly. 


CHAPTER 2 
Software 


2.1 Software Overview 

The original computer program was written in Microsoft Quick C (almost 
ANSI C compatible) to run on a PC-386; C was chosen for several reasons, the 
main ones being, - . — 

1. it has a relatively fast compiled code execution, 

2. it is the most commonly used language amongst the students 

3. and faculty working in the laboratory 

4. it is a flexible language capable of implementing complex data structures. 

The PC-386 has now been replaced by a set of four micro-processors, the Inmos 
Transputers, interconnected by high speed serial data links and a host computer (a 
PC-286) for file handling and the user interface in order to greatly increase compu- 
tation power and to facilitate further expansion if so desired; the Transputers and 
their links constitute the computer hardware network (see figure 2.1). 

Optimum performance is obtained if software is written in the language Occam 
which was developed by Inmos for the Transputers and therefore the speed critical 
control and data 1/ 0 interface algorithms were converted to two separate processes 
using the Inmos Parallel Occam compiler while the rest of the original program (i.e.. 
the user interface and file handling section) was recast as a single process using 
the Inmos Parallel ANSI C compiler. The processes communicate with each other 
and the host through asynchronous data channels; the processes and their channel 
configuration constitute the software network and this is distributed amongst the 
hardware network as shown in figure 2.1. 
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The user and file interface process is called flexl.c and it provides data input 
and storage and user communication via the host computer screen and keyboard. 
Data files are stored in the Mathworks MATLAB .mat format to allow an easy 
interface with the tools of the MATLAB mathematical analysis software on either 
the host PC-286 or a SUN workstation network to which it is connected. Figure 2.2 
shows a high function level schematic for flexl.c which emphasizes the .mat file 
interface and the initialization and/or resetting sequences which the testbed is put 
through prior to each experiment. 

The process ctalg.occ implements the several realtime control algorithms used 
during experimentation. These are written as general discrete time state space sys- 
tems and/or discrete time nonlinear operators whose parameters are configured prior 
to am experiment from data stored in a .mat file. Input consists of a combination of 
sensor feedback from the testbed and pre-defined data vectors while outputs are the 
actuator control signals (only one is needed to control the single link flexible beam). 
The two control algorithms used in this report are, 

1. Open Loop Control which passes input data from a .mat file directly to the 


actuator channels. 
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Figure 2.2: Software Flow Diagram #1 (High level) 
2. Linear Observer and Feedforward Compensator. 


.2 Program Description 

The main purpose of the software is to provide a set of easily configurable 
constant time interval sampled data controllers for the flexible beam apparatus 
(i.e., the “plant”), including a linear observer/feedforward compensator, with the 
following features (see figure 2.2), 

1. Each controller has a fixed structure written in software and selectable pa- 
rameters which are loaded from a .mat configuration file. Note that there is a 
compromise to be made between the generality of the controller and the speed 
of execution, e.g., the linear observer/ feedforward compensator controller may 
be written as a more general linear controller however matrix computations 
increase. 

2. Controller input and feedforward trajectories are stored in a .mat input file 
as data vectors (i.e. 1 x n, or n, x 1 arrays where n, is the number of data 
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samples taken during the experiment), 

3. Controller output (i.e., plant input) and plant feedback to the controller vec- 
tors are stored in a .mat output file after the experiment. 

4. Plant input and output signal steady state offsets may be automatically mea- 
sured prior to each experiment, stored in a .mat file, and compensated for 
during the experiment. 

5. Both the sampling time interval t s and the number of samples taken during 
an experiment are specified in the input file. 

2.3 Resetting Signal Offsets 

Each analog component of the testbed hardware is subject to nominal signal 
changes between experiments due to temperature induced drifts, aging of compo- 
nents, and distortions caused by the experiments themselves (e.g., a large strain on 
the flexible beam can slightly bend it so that the strain gage steady state values mea- 
surably change). In order to minimize the effects of such errors the software includes 
an option to measures these offsets prior to the experiment and the I/O procedure 
then automatically compensates (by subtracting from A/D values and adding to 
D/A values). Furthermore the new offsets are stored in a .mat file (mfl.cfig.mat) 
and are used during subsequent experiments until they are measured again; in prac- 
tise the offsets are usually updated prior to each experiment. 

Strain gage signal offsets are measured by averaging a series of readings taken 
over five seconds while the beam is at rest (i.e., the input torque is at zero) while 
the potentiometer offset is always set to zero for experiments where the absolute 
zero angle is unimportant. 

The offset related to the motor amplifier circuit is more problematic due to 
the effect of hub shaft friction. In the absence of friction the unique amplitude of 


14 


the signal applied to the power amplifier would equal the offset when the motor 
is at rest and hence be easily measured. However in actuality the beam will come 
to rest if the applied motor field torque is less than the friction dynamic threshold 
and remain there as long as the static threshold in not exceeded and hence zero 
applied torque condition is not readily detected. One solution is to apply a dither 
signal directly to the motor amplifier and add a proportional negative feedback loop; 
assuming a sufficiently large dither signal and convergence of the operating point to 
a periodic cycle the amplifier circuit offset will correspond to the average applied 
torque signal. In practice a 10 Hz signal is applied and after a pause of 5 seconds 
the average over 5 seconds of the torque signed is calculated. 


CHAPTER 3 
Beam Modeling 


3.1 Beam Analytical Modeling 

3.1.1 Linear Differential Equations 

The essential model components of the testbed are the beam itself driven by 
a torque at the hub, and the effective hub inertia. Figure 3.1 shows the coordinate 
frames and geometric notation used in the following equations. The correspond- 
ing non-linear beam model was derived in [1] using Hamilton’s Principle with the 
following initial assumptions, 

1. Hooke’s Linear Stress-Strain Law holds. 

2. The beam elongation is negligible. 

3. The beam deformation is small. 

4. The hub and bending velocities are small. 


and a set of linearized partial differential equations was subsequently obtained, 


8 A v 8 2 v 
EIa 8x< +Px Ji* 

T-I H j> + EI a v"{0,t) = 0 


(3.1) 

(3.2) 


with the boundary conditions 

u(0, t) = 0,/; u'(0,f) = 0, v"(0,i) = 0, v"\L,t) 
E = beam material Youngs Modulus 
r = the input motor torque, 
p x = the beam mass per unit length. 


= 0 , 
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Figure 3.1: Beam Model Coordinate Frames 


Ik = hub inertia, 

L = the beam length, 

V = U) + X<t>, 

4> = rigid body mode angle, 

Oh = the hub angle, 

u>(x) = beam displacement at x meters along the nominal axis. 


3.1.2 Normal Mode Form 


An eigen-analysis of these equations is performed in order to establish the 
infinite set of modal frequencies {u>,} at which the beam exhibits resonance (damping 
is not included in the model at this stage) and corresponding mode shapes {^\(x)}, 
including one rigid body mode (wo = 0, tf> o(x ) = a straight line) where, 


^i(x) 

Ai 

Bi 

Di 


Ai sin(fcix) + Bi sinh(fc,x) + cos(fc,x) + D n cosh(^x) 
c ch + $ sh + 1 


2 p x 

hkf 

-C, 


•s ch + csh 


-Ci 


C, - Ai 


(3.3) 

(3.4) 

(3.5) 

(3.6) 
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Computer + 0/A & A/0 Hardware 


1. Linear model 

2. Low frequencies 


f Q {x) = x and }f{t,t) = t. 
T am (s) ~ 1 and T ssc (s ) ~ 1333 . 


Figure 3.2: Testbed Physical System 
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Figure 3.3: First Four Normal Mode Shapes 


where s = sin(fc,Z-), c = cosfaL), sh = sinh(fc,£) and ch = cosh (&,£) 


C, may be chosen arbitrarily and may be used to scale the units of {t/>,(x)} 
(dimensions meters). {A:,} are the set of solutions to 

—k?{l+cch) + {sch-csh) = 0 (3.7) 

Px 

and the corresponding oscillation frequencies are given by, 

= kjJ^- (3.8) 

Table 3.1 lists the first nine such frequencies for the flexible beam used (see 
appendix 5.3 for the material dimensions and constants) plus a zero which represents 
the zero k solution to equation ( 3.7) and the rigid body mode; figure 3.3 shows the 
corresponding beam mode shapes described by equation ( 3.3) and scaled so that 
each a n is unity. 

Note that the mode number is the same as the number of oscillation zeros that 


exist along the beam. 
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Table 3.1: Beam Normal Mode Parameters 



1. * = mode number (i.e., number of mode shape nodes — 1) 

2. fi = calculated mode natural frequency = ^ 

3. fi = measured mode natural frequency = 

4. r, = measured decay time constant = 

5. = measured damping ratio 


Next the linearized partied differential equation solutions axe obtained in terms 
of these modal quantities to yield v(x,t), the linear approximation of the radial 
displacement of the beam at a distance x along the beam from the hub at time t, 



oo 

v{x,t) = £$,•(<) V>i(x) 

i— n 

(3.9) 


- 2 «(0) r 
$. + qi = .r 

(3.10) 


<HPx 



where qft) = the i tk modal amplitude function , 


* ; 

a, = < > +(— )#(0)$(0), 

Px 

(3.11) 


<f(x),g(x)> = i2[0, L], the innerproduct of f(x) with ^(x), 



The ith mode is thus governed by a second order linear differential equation 
( 3.12) which is an undamped forced harmonic oscillator. The pole damping ratio 
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(£,-) for each flexible mode is now estimated experimentally and inserted to produce, 


qi + 2C,a'i?i = 


W) 

0>ipx 


(3.12) 


3.1.3 Finite State Space Normal Mode Approximation 


The A and B matrices of a finite state space representation n modes can now 
be easily derived from equation ( 3.12) if we truncate the model after the n th mode 


x = 


A m x + B t t where, 


(3.13) 



0 

1 

0 

o ... 

0 

0 


0 


0 

-do 

0 

o ... 

0 

0 


*n(0) 
ap i 


0 

0 

0 

i ... 

0 

0 


0 

A m — 

0 

0 - 

-ul 

— d\ 

0 

0 

, Br = 

*J(0) 

txp x 


0 

0 

0 

0 ... 

0 

l 


0 


0 

0 

0 

0 


-dn 

1 T 


^(Q) 
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qo 

<io qi 

q i 

••• qn - 1 

qn- 1 




(3.14) 


(3.15) 

Currently used are a hub angle sensor (a potentiometer) and four strain gages 
mounted along the beam and in each case the signals produced can be approximated 
by a linear combination of the states x(i ) to produce the output C matrix, and the 
D matrix is null and the corresponding sensor velocities are obtained by shifting the 
appropriate C matrix rows one place to the right in order to meet the q n states. 

3.2 Hub Actuator 

With reference to figure 3.2 the switch mode power amplifier and DC motor 
dynamics may be approximated by the following linear transfer function relation- 
ship, 


r(s) = k T k am T am (s)V m (s) 


(3.16) 
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where r am (s) = 
and k T = 


Tim 
“I" Tim 
T ( s ) 

- V , the motor torque constant. 

AnOO 


(3.17) 

(3.18) 


Here r am is the dominant time constant of the power amplifier and DC motor 
field winding circuit, T am (s) is the associated low pass transfer function due to 
resistance and inductance (see chapter 1) and k am is the low frequency switch mode 
power amplifie/DC motor voltage to motor circuit current ( I m ) gain. Since r am 
and is estimated to be 2 milliseconds approximately; it is neglected (i.e., assumed 
infinite) for the purposes of this report as the highest modal frequency considered 
is 58 Hz, and hence equation ( 3.17) may be approximated by 


r(s) — k r k am V m ( K s') . 


(3.19) 


3.3 Hub Angle Sensor 

Equation ( 3.13) describes the torque to modal state transfer function and we 
now relate the modal state to the hub angle <j> ; the truncated form of equation ( 3.9) 
implies 


*W(0,i) = EU*(Otf(0) (3.20) 


where i/>'(0) is the derivative of the i th mode shape at the origin (i.e., the hub 
shaft) and is found by differentiating equation ( 3.3) . Equation 3.20 can be recast 
in state space form as, 


<t> 

c V 


C+x 

o v-i(o) o ... v;to) o] 


(3.21) 

(3.22) 


Note that C * can be related to B r (see equation ( 3.14) as follows: 


CV = 


(3.23) 
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where P 4 is a diagonal strictly positive definite matrix and since the only 
system pole that is not exponentially stable (i.e., the rigid body mode pole at the 
origin) is decoupled from the hub velocity output the Positive Realness Theorem 
guarantees that the torque to hub angle transfer function is strictly positive real, 
(see chapter 4). 

3.3.1 Angle to Voltage Conversion 

A single turn potentiometer is used to convert the hub angle <f> to a voltage 
signal: with reference to figure 3.2 

V a = <i>.k v and k p = ^ (3-24) 

where V s is the potentiometer supply voltage, and k p is the conversion constant. 


3.4 Beam Strain Sensors 

The longitudinal strain at a specific point on the beam surface is closely related 
to the curvature w"(x, t) of the beam at that point and we can relate the latter to 
the modal states in truncated form as follows, 

w"(x,t) -»'(*,«) = £ «.(!)#'« (3.25) 

1=1 

where ^'(x) are given by equation 3.3 and therefore in state space form, 


w 


"(*,*) = C,(x)x 


C,(x) = 


0 0 V>"(x) 0 ... i>"(x) 0 


(3.26) 

(3.27) 


Note that i/> 0 (x) may be calculated explicitly and is found to be 

1 r \ r/L 3 , ^ H \ 
tl> 0 (x) = Va(— + —)x 

o p x 


(3.28) 


and hence tPq(x) is zero for all x. 



23 


3.4.1 Strain Gage Geometry 

The strain gages are mounted in pairs opposite to each other on either side of 
the beam so that together they can measure the extent of bending at a particular 
distance along the beam, half way up the side, (see figure 1.1). Only one strain 
gage is strictly necessary at each point of interest, however considerable noise and 
disturbance immunity is achieved by using them in pairs as one arm of a Wheatstone 
bridge configuration (see figure 3.6). In fact four strain gages completing the bridge 
circuit are even better, two mounted on each side of the beam close together, but it 
is not clear that the performance improvement over two is significant and worth the 
extra cost and construction. Figure 3.4 shows a cross section of the beam with a pair 
of strain gages on either side (looking down from above the beam). The important 
dimensions are, 


dbeam — the width of the beam, 
d,ub = strain gage substrate thickness, 
dgiue = estimated glue thickness, 
dfiim = thickness of strain gage film, 


and di /2 


— dbeam /2 “ 1 " dg[ uc * 4 “ d 3U b 4 * djUmj 2 

« 7.938 x 10 -4 + 5 x 10 -7 = 8.438 x 10 -4 


(3.29) 

(3.30) 


is the distance between the beam neutral axis and the center of the strain gage 
film half the effective beam width with respect to the strain gages. 

A strain gage, of the type used, is a pattern of resistive material whose resis- 
tance increases when stretched and when the beam bends in a clockwise direction 
(cw) (looking down from above) the cw+ strain gage is elongated and so its length 
undergoes positive strain (hence the name cw+), while at the same time, the cw- 
strain gage is forced to contract, implying negative strain. The strain gages are very 
linear devices over a wide range of strains and the altered resistances increase V sg 
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Figure 3.4: Beam and Strain Gage Cross Section 



Figure 3.5: Beam Curvature Geometry 


in figure 3.6 which is amplified by the strain gage signal conditioner to produce V so 
while the opposite occurs when the beam curves the other way. 

Note that the strain gages directly measure strain in the plane of the strain 
gage, not on the surface of the beam or elsewhere, however the strain in any plane 
parallel to the nominal axis plane is linearly proportional to its distance from the 
nominal axis plane, and hence conversion between different strains is achieved by 
scaling. 

The quantities of interest in this project are the modal state amplitudes and 
we can relate these to the strains by first constructing figure 3.5, a cross sectioned 
diagram in the x-y plane of the rotating beam frame. Let x\ an arbitrary point 
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along the beam neutral axis when at rest (ie. the x axis) with iw x (= u;(x x )) as the 
nominal axis displacement associated with the linear model; n x = (xi,u;i). Also 
choose x 2 such that 


x 2 = Xi + Ax, 

tv 2 = w(x 2 ) ;n 2 = (x 2 ,ty 2 ). 

It is assumed that Ax is chosen sufficiently small such that the segment of the 
n.a. C n , ( = [ n x , n 2 ] ), is continuously convex w.r.t. a particular side of the 
beam. Let A n be the corresponding nominal axis curve section, length /„ (i.e., A n is 
a circular arc centered at p and connecting n x to the transverse plane line through 
n 2 ). At n x and n2 extend the beam transverse lines until they meet at some point 
p and let R, be the length of [ n, , p ] ; 0, is the angle between the x axis and the 
transverse lines. The transverse plane and outer surface intersections axe labeled s x 
and s 2 and the connecting beam curve is called C a with a corresponding nominal 
axis curve section A s . 

Let 


9 ti 

= 

&i — 90 degrees , * = 1,2 (then w\ = tan(^ti)) 

(3.31) 

Ad 

A 

0 2 — 0\ — 6t 2 ~ dn, 

(3.32) 

1 cn 

= 

length(C„), /an = length(A„) 

(3.33) 

L c , 

= 

length(Cj) = /cn + Al C s 

(3.34) 

L a5 

= 

length(A,) = l an + A l aa 

(3.35) 


where A and A l an are the n.a. curve and arc length extensions due to strain 
deformation respectively. 
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3.4.2 Curvature-Strain Linear Relationship 


The following limits are assumed 


Now 


,. Cn 1 1 . 

lim . — 1 5 lim . — 

Ax-0 A n Ax— >0 A, 

lim — — 1 , lim -p— ■ = 1 , 

Ax— *0 l an Ax-0 L a!> 

.. A/ OJ .. //» , 

lim — t- = 1 , lim - — = cos( 0 n )- 

Ax— ►O A/ c . Ax— ^0 l C3 


Rl A6 = la n* (i?i + ^p)A0 = / on + A/ OJ , 


dfteom 


A0 = A/ a 


(3.36) 

(3.37) 

(3.38) 


(3.39) 

(3.40) 


i-i# ^ . 4 i 


and therefore 
*. 

W 2 — w[ = sm(A0)/(cos(O n )cos($t2)) 


n 1 - (**>1 - w i) r 

;uuk.- so , = lim 7 = hm 

1 Ax— 0 Ax 


sin(Ad) 


1 


lim£ = 


Ax— 0 cos(^a) cos(0 (2 ) Ax 

A la 


COS 2 (#ti) Ax— 0 Ax 


lim 


1 


lim 


AL 


COS 3 (0«) Ax-0 /, 


COS 2 (0 tl ) Ar— 0 Ax 
1 

lim 


1 


lim 


AL 


AL 


d itam cOS 3 (0 tl ) Ax— 0 / c 
ei 


COS 3 (0 t i) 


cos 3 (0 a ) 


(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 


where ei is the beam convex surface engineering strain (and a convex surface 
implies that ej is positive). Lastly 


A .. A lea 

ei = hm — — 

'«-* 0 la , 


and hence e,i = ei-7 3 — — w"( cos 3 (^n)) 

“1/2 * 

=► Ulj ~ ( ; )C 4 J 

«l/2 


(3.46) 

(3.47) 

(3.48) 


for small deflection angles 9 n where e, is the strain experienced by the corre- 
sponding strain gage, (the strain in the opposite strain gage is equal to — e,). The 
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term — r 2 — is due to the fact that the active part of the strain gages are raised from 
<* 1/2 

the beam surface and hence the true distance from the center of the beam is larger 
than halve the actual beam width. 


3.4.3 Curvature-Strain Non-linear Relationship 

Equation 3.47 is nonlinear and may be expanded as follows 

w\ = tan(0 a ) ^ cos 2 (0 t i)(wi) 2 = (1 — cos 2 (0 a )) 
=» cos 2 (0 n ) = (1 + (twi) 2 ) -1 

=> cos 3 (0 n ) = (1 + (^i) 2 ) -2 ^ 3 

d b 

=> = 


a. < 


( i + K) 2 ) 2/3 ' 

Assuming that {w[) 2 « 1, a binomial expansion implies that 

dfceom 


- 


L <(l-3/2 K) 2 + 0(K) 4 )) 


and since w[ may be expressed in terms of the modal states 


u>i(0 = v'(x u t) - <f>(t) = JZ ~ ^i(O)] 

i=0 


(3.49) 

(3.50) 

(3.51) 

(3.52) 


(3.53) 


(3.54) 


equation ( 3.48) may be modified to included a power series approximation to 
the nonlinear relationship up to any desired order : 


w. 


iT-h 


•e.i- 


(3.55) 


l </i/ 2 '(l-3/2W) J + 0(Wn) 

For the discrete model a previous estimate of the plant may be used in order 
to calculate w[ and this then be subsequently used in equation ( 3.47) for the 
current estimate of w ". This approximate non-linear compensation is not used in 
simulations or experiment in this Thesis. 
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3.4.4 Strain to Voltage Conversion 

As mentioned above two strain gages are used opposite each other at every 
measurement point along the beam in a wheatstone bridge circuit, (see figure 3.6). 
The 10V DC supply voltage ( V , p ) is provided by the Analog Devices 3B18 signed 
conditioner and the bridge output ( V ag ) feeds into this same circuit which amplifies, 
buffers and low pass filters the signal (see appendix 5.3). If Rcw+ and R^u- are the 
resistances of the strain gages cw+ and cw- respectfully (see above) then the basic 
strain gage equations are 

= Ke cw+ , = Ke cw . = -A'e w+ (3.56) 

1LCW+ Hew— 

where K = some constant dependant upon the strain gage selected and A R cw+ 
is the resistance change in strain gage cw+ due to strain €<*„+. 

In figure 3.6 Rn = Rcw+ nominal = R^- nominal, and R + , R- are chosen 
equal to Rn also so that when the beam is at rest V, is zero. Then 

V, = and V 2 = (3.57) 

=5- V„ = V,-V, = (3-58) 

Combining equations ( 3.48) and ( 3.58) gives (for each strain gage pair) 

V Hl = (3.59) 

where k s = — ^ — -di/ 2 - (3.60) 

Note that in figure 3.2 V, g i is called W^. K = 2.08, V ap = 10 and d ly r 2 as 

given in equation ( 3.30) implies that k, = 8.775 x 10 -3 . 


3.4.5 Strain Gage Signal Conditioner 

The strain gage signal conditioner contains a second order low pass filter given 
by 

A 

T - k — 

s 2 + 2Qssc^ssc + ^isc 


(3.61) 
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(3.62) 

Currently k 3ac = 1333, C«c = 1 and ui„ c = 100.2.x however these dynamics are 
not included in the linear model used in this Thesis, rather only the low frequency 
aproximation of T MC (s) ~ k slc is used. 

3.5 Computer Input and Output 

Each A/D converter may be considered as a constant gain of k a d (ad units )/ V 
followed by an amplitude quantization operator /q(.), while a D/A converter is the 
is the same operator followed by a constant gain of k^a. which is currently configured 
so that kda = -A-, i.e., 

K ad 

y ad (t) = f Q (k ad x v (t)) (3.63) 

and y v (t) = k da fq{xd a {t)) (3.64) 

where x„(t) and y v (t) are voltage signals, and x^t) and y a d(t) are amplitude 
quantized signals. Jq(x) rounds x to the nearest integer but this this quantization 
effect is ignored in the linear model and so 


J/oci(0 — k a dXv(t)i 


( 3 . 65 ) 
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Vv{t) = k d *Xda{t) 


(3.66) 


are assumed. 


3.6 Complete Testbed State Space Model 


With reference to figure 3.2 the separate system transfer functions may now be 
cascaded to produce the overall model referenced to the computer control algorithm 
as follows. 


x 

2nx2n 


A 

g2nxl 


Ax + Bu, y = 
A m 

B 7* ^7 Jc a 

kpC -[p 

k v Cl 


Cx 


>2mx2n 


— kad 


kssc k* C$q 

L L nv 

H'SSCA' 3 '^' sQ 


where x = 


u = 


(3.67) 

(3.68) 

(3.69) 


(3.70) 


and assuming y = 


C v 

^ row 


kssckgC s4 

k k C v 

?0 4o <?1 4l ••• ?n- 1 1 ] , (3.71) 

control effort in AD units (3.72) 

j/Au6 yLb y*o I i v ,o - y * 3 y v s3 ] in AD unit * 3 - 73 ) 

0 C row (l,...,2m-2) ] (3.74) 


i.e., C T ow shifted once to the right and with a zero inserted. 

Figure ( 3.7) shows the bode plots for several input/output transfer functions 
when the first five flexible modes are included; it is emphasized that the above 
model, and therefore these plots, include damping information that was obtained 
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experimentally, however the form of the plots is invariant wit respect to large es- 
timated damping ratio errors. Note that the hub angle velocity transfer function 
bode phase plot is bounded by ±90 degrees as expected because it is strictly positive 
real, however so is the phase plot for the first strain gage response even though the 
sensor is located a distance of 6.5mm away from the hub along the beam and this is 
explained in the next section. All the other strain gage velocity transfer functions 
are seen to be not positive real. 

3.7 Model Structural Features 

Figure ( 3.8) includes a schematic representation of the plant model which 
helps to illustrate the following points. 

1. The single input torque simultaneously excites a low pass filter and several 
second order filters, while each sensor is a linear combination of the filter 
states. 

2. The rigid body mode transfer function, torque to qo is a passive low pass filter, 
G 6 (s) = s/(s + r). 

3. Each flexible mode transfer function, torque to q n , n > 0, is a passive bandpass 
filter, G;(s) = 2 £;u>;s/(s 2 + 2£,u;;s + uif). 

4. The torque to hub angle velocity transfer function is a positive linear combi- 
nation of the G^(s) an< ^ h ence it is minimum phase. 

5. In the complete linear model, each torque to strain velocity transfer function 
is not strictly positive real and hence there exist negative C s (x) elements 
(assuming x > 0). In general the further the strain gage is located from the 
hub the “sooner” the first negative C,(x) element occurs for a given number 
of modeled modes or, alternatively, for a given non— zero distance from the 
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Figure 3.7: Finite Linear Model Bode Plots 
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Plant 



Figure 3.8: Normal Mode Linear Model Structure 


origin, a negative element will occur in C,(x) for a sufficiently high order 
model. However, if the strain gage is close enough to the hub, C s (x) may 
contain only positive entries for any given number of beam modes and this is 
the case for the first strain gage (#0), which is only 6.5mm from the hub, and 
at least the first five flexible beam modes. 

3.8 Model Identification 

A modal analysis was carried out on the testbed in order to identify the three 
sets of finite linear normal mode parameters, i.e., 

1. the first n modal oscillation frequencies (/ n ), 

2. the corresponding pole damping ratios (£„), 

3. the corresponding input/output gains for each sensor 

and the results for /, and £,• are shown in table 3.1 for 0 < i < 5. Modal analysis 
was used since the test signals may be selected to target individual parameters of 
interest, e.g., a sinusoid at a specific flexible mode frequency produces a response 
strongly dominated by the corresponding flexible mode and the input/output gains 
may now be directly determined. The identified model parameters are considered 
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more accurate than the analytically derived ones from the beam and hub assembly 
physical parameters since some of the latter are only known approximately however 
the analytically derived input/output gains gains are used in the final model (i.e., 
the model used for controller design and analysis in the next chapters) since the 
identified input /output gain data is not up to date. 

3.8.1 Oscillation Frequencies 

The first step was to apply a narrowband test signal such as a PRBN function 
to the beam and obtain an FFT of any sensor response (although the strain gages 
are preferred as the oscillation peaks are well defined) and hence a set of bode plots. 
Since each mode is very poorly damped the shape of the sinusoidal response of a 
sensor close to an oscillation frequency is strongly determined by the corresponding 
modal transfer function and hence the experimentally obtained bode magnitude plot 
peak locations determine the oscillation frequencies (and hence the associated pole 
magnitudes). A second refined bode magnitude plot was then obtained by applying 
individual sinusoids at frequencies close to the estimated pole magnitudes and then 
the latter were re-evaluated. 

3.8.2 Damping coefficients 

The rigid body mode non-zero pole was measured by applying a rectangular 
input torque and curve fitting the hub angle response in order to determine the rate 
of exponential increase (a rectangle rather than a step signal was used since the 
latter would of course cause the beam to wind out of its operating envelope). This 
mode does not affect the strain gage responses and hence the corresponding strain 
gage data cannot be used in order to estimate rigid body mode parameters. 

In order to measure the i t>l flexible mode pole damping ratio (Ci) a sinusoid was 
applied at the corresponding center frequency (u>i) until a steady state situation was 
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reached and then abruptly removed at time fj; the envelope of each sensor response 
after <1 exponentially decays with a time constant f,- given by 


Ti = 


1 

C.A- 


(3.75) 


and once t,- is measured by curve fitting this envelope may be calculated 
since u>i is known. 


3.8.3 Modal Gains 

This last set of data also yields the input/output gains 6,c,j for each sensor j 
at each oscillation frequency u>,- from which the state space model B and C matrices 
may be formed. 

3.8.4 Identification Improvements 

The above procedures may be improved in order to more accurately identify 
system modal parameters and some ideas are as follows. 

1. A discrete fourier transform is an estimate of the continuous fourier trans- 
form at specific discrete frequencies that are determined by the data sampling rate 
and the number of samples. This suggests that a better narrowband signal than 
PRBN would be one composed of equal amounts of sinusoids at the discrete fourier 
series frequencies with a random constant phase shift applied to each sinusoid in 
order to reduce coherent peaks in the time domain signal: unlike PRBN this sig- 
nal does not waste energy on frequencies other than those assumed by the discrete 
fourier transform (and therefore those of interest) and provides equal weights for 
all frequencies. Furthermore, the amplitude weights may be accurately altered if 
desired, perhaps to increase the measurement signal/noise ratio at frequencies close 
to plant transmission zeros and any number of sinuoids may be omitted in order to 
precisely tailor the applied signal bandwidth. This test signal was used for the Bode 
plots shown in figure 5.11 which are much smoother than those obtained with PRBN 
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test signals. It may also be possible to specify the individual sinusoid constant phase 
shifts rather then set them randomly in order to minimize the maximum test signal 
peak in th etime domain so as to further avoid nonlinear response effects associated 
with large input signals. 

2. Estimate the coulomb friction and include corresponding nonlinear com- 
pensation during the identification process. This is particularly important since 
test data show significant nonlinear effects when small signals are applied to the 
beam and these are assumed to be due to coulomb friction. Furthermore, very large 
applied signals also result in nonlinear effects due to nonlinear damping, stiffening 
and gravity /torsion effects which may be more easily avoided if smaller signals are 
allowed. 

3. Calculating complex pole damping ratios by analysing sinusiod decay en- 
velopes can be difficult due to nonlinear effects that add unwanted frequencies to 
the data and one possibility is to filter the data with a linear second order band- 
pass filter centered at the same modal center frequency first and then the modal 
damping ratio may be derived from the measured value and that of the bandpass 
filter: the center frequency for such a filter is known exactly since it corresponds to 
the applied sinusoid frequency. Another option is to filter the data in a nonlinear 
manner by applying a window to the corresponding discrete fourier series to select 
those frequencies close to the modal center frequency and then transform back to 
the time domain in order to analyse the decay envelope : the decaying sinusoid 
fourier transform is dominated by frequencies close to its center frequency. 

4. The rigid body mode non-zero pole was identified with a rectangular test 
signal which excites the flexible modes also and though this was taken into account 
a better result might be obtained by using a smoother signal, perhaps a rectangular 
pulse with parabolic ramping or filtered with a sharp roll off low pass filter. 


CHAPTER 4 
Control Design 


4.1 Overview 

The linear modal flexible beam model developed in chapter three (see equation 
( 3.67)) possesses a simple decoupled structure , i.e., a collection of second order 
systems (a single rigid body mode and the set of flexible modes) driven by a common 
input torqe. The rigid body mode largely determines the step response time while 
the flexible modes govern the oscillatory behavior and applying feedback control can 
decrease the step response time (i.e., move the rigid body mode real poles away from 
the origin in the Laplace s domain) and increase the flexible mode damping ratios 
(i.e., swing the flexible mode complex pair poles towards the negative real axis in 
the Laplace s domain). Applying any amount of hub angle proportional feedback 
moves the pole at the origin and the new closed loop system becomes completely 
controllable and observable, hence full state feedback may be used to arbitrarily 
place the closed loop poles. However in reality noise, non-linearities and modeling 
errors limit the extent to which feedback may be applied and the problem now is 
how to select a robust feedback controller. 

A useful design tool is the Passivity Theorem (see reference [4]) which guar- 
antees that strictly positive real feedback compensation of a positive real (m inputs 
and m outputs) plant will produce an asymptotically stable closed loop positive real 
system even in the presence of modeling errors, so long as there are no pole-zero 
cancellations, and furthermore such a compensator will have an infinite gain margin 
when applied to the plant model. In the case of the flexible beam a natural choice 
for a positive real (with respect to the input torque) output is the time derivative 
of the hub angle which may be obtained with an appropriate sensor, direct filtering 
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of the hub angle, or reconstructed from an estimation of the plant states (i.e. using 
a closed loop observer). The Passivity theorem may now be applied to show that 
hub angle proportional plus derivative (PD) feedback produces an asymptotically 
stable closed loop positive real system with an infinite compensator gain margin. 
This is the basis of the first controller designed and root locus analysis is used in 
a systematic fashion in order to select the actual proportional and derivative gains 
as a compromise between rigid body mode performance, flexible mode performance 
and pole location sensitivity to gain variations. 

However the second and higher order flexible modes are only weakly observable 
through the hub angle sensor and in order to substantially increase the damping 
ratios of these modes requires unacceptably high PD gains. A solution to this 
problem is to synthesize a new positive real output signal, either directly from 
several plant sensors (collocated and/or non-collocated) or via an estimate of the 
PD compensated plant states, that is strongly coupled to the modes of interest 
and then apply a second positive real compensator (see figure 4.5). The Positive 
Realness Theorem provides a means for constructing such an output from the PD 
compensated plant states and a second controller is designed this way that improves 
the damping ratios on the second and third flexible modes using static feedback; a 
third controller tested uses a positive real filter in the feedback loop in order decouple 
the dynamics of higher order uncompensated modes. 

Each of these controllers is designed using a four mode plant model (eight 
states) and an observer to provide an estimate of the plant states from which positive 
real outputs are generated. The effect of some unmodeled higher order modes on 
the positive realness of synthesized plant outputs is examined (through simulations) 
with different closed loop observer types, namely, full order pole placement, reduced 
order Luenberger pole placement and LQE designs and the observer which produces 
an output least sensitive to the modeling error is selected: open loop observers are 
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not tested. 

When designing and testing compensators a mixed analytical/identified model 
was used since the identified model parameters are considered more accurate than 
the analytically ones derived from the beam and hub assembly physical parameters 
(see chapter 3). 

4.2 Positive Realness 

With reference to [2], assume a minimal (controllable and observable) expo- 
nentially stable system S with time domain representation 

x(t) = Ax(t) + Bu(t), x(0) = x 0 

y(t) = Cx(t) + Du(t ) (4.1) 

and frequency domain representation 

y(s) = (D + C{sl - A) -1 B)tz(s) = T(s) (4.2) 

where x € u, y € and cr m i n (B) > 0. 

Theorem la : 

(A, B, C, D) is Strictly Positive Real if there exists P > 0 , P , L € 3? nXn , 
fimin(L) = € > 0 , Q € 3R mxn , and W € that satisfy the Lur’e equations 

A t P + PA = -Q T Q - L (4.3) 

B t P-C = W t Q (4.4) 

W T W = D + D T . (4.5) 

For the flexible beam model (3.67) D is null and so we may take W and Q to 

be null in which case the above three equations reduce to 

A t P + PA = —L 

= b t p. 


c 


(4.6) 

(4.7) 
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This theorem offers a method of computing an output matrix C in terms of 
L such that the corresponding system is positive real as follows : equation 4.6 is a 
Lyapunov matrix equations and since fi max (A) < 0 (e.g., S is a minimal exponen- 
tially stable system) there exists a unique solution P for any given L (see [5], page 
574) which may be expressed as 

P = J\ ATt Le At (4.8) 

Furthermore, 

L > 0 and L* = L =>• P > 0 and P* = P. (4.9) 

Therefore B T P(sI — A)~ l B will be a positive real transfer function if L is 
restricted as in 4.9. In practice equation 4.8 is not suitable for computer computation 
but other more efficient methods exist and are readily available (e.g., the LYAP 
function in MATLAB). 

Physically, L determines the convergence rate of the Lyapunov function V (x) = 
x m Px w.r.t. the open loop plant since 

^F(z(f)) = —x*(t)Lx(t) (4.10) 

along any solution trajectory of x = Ax. 

Now, if D is null, the equivalent frequency domain conditions for positive 
realness are: 

Theorem 2b : T(s) is Strictly Positive Real if for all w € 3? 

T(juj) + T*(ju>) > 0 (4-11) 

and Jirn + T*(ju)) > 0 (4.12) 

Hence a Nyquist plot for a positive real transfer function will lie in the closed 
right half complex plane while a bode phase plot will be bounded by ±90 deg. All 
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Figure 4.1: SISO System 


minimum phase first order filters, such as Po{s), are positive real and so are second 
order bandpass transfer functions, such as P n {s)- Furthermore, 

Pi(s) positive real and C; positive constants for all i € {0, 1, ...(n — 1)} 

=*■ Re ^ Pi(s)Ci = Ci[Re P,(s)] > 0 => ^( 5 )C* 13 positive real, (4.13) 

i=0 «ssO »=0 

i.e., positive linear combinations preserve positive realness. Since the hub angle 
velocity output of the flexible beam linear model is a positive linear combination of 
the modal velocity states (themselves positive real outputs, in fact strictly positive 
real outputs, with respect to the input torque) the torque to hub angle velocity 
transfer function is positive real. A similar result is true for strict positive realness. 

4.3 The Passivity Theorem 

With reference to figure 4.1 

Theorem 3: 

If P(s) is positive real and C(s) is strictly positive real then 
T(s) = P(s){sl - C(s)P(s))~ l is BIBO stable . 

Furthermore, if there are no unstable pole zero cancellations then 
T(s) is internally asymptotically stable. 

This is a particular case of the Passivity Theorem and can be applied in the 
case of hub angle PD negative feedback for the flexible beam to guarantee closed 



loop asymptotic stability assuming a non-zero derivative gain and more generally 
any collocated actuator and velocity sensor pair will produce a positive real transfer 
function for a mechanical structure that includes no energy sources. 

Note that the Passivity Theorem is sufficient but not necessary; less restrictive 
conditions based upon the open loop plant may be developed from the Nyquist 
stability criterion which is itself a necessary condition. 

4.4 Hub Angle PD Compensator 

The hub angle forward path PD compensator is designed in two stages; firstly 
only the rigid body mode of the plant is modeled and a two column matrix of 
increasing PD gains is obtained analytically corresponding to increasing closed loop 
step response performance. Secondly the modeled flexible modes axe included and 
a particular pair of PD gains is selected based on the characteristic equation root 
locus plot. While pole location sensitivities w.r.t. compensator parameters are also 
considered no attempt is made to analytically or numerically optimize a performance 
and/or robustness index. 

4.4.1 Closed Loop Stability 

Figure 4.2 (a) shows the flexible beam plant (torque to hub angle) with a 
PD feedback compensator. By transferring the plant integrator to the compensator 
block as in part (b) of the figure the modified plant becomes strictly positive real 
while the compensator becomes positive real (and strictly positive real for d > 0) 
and hence the Passivity Theorem guarantees asymptotic stability of the the closed 
loop system. 
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Figure 4.2: PD Compensation Systems 


4.4.2 Rigid Body Mode Plant 


Consider the closed loop transfer function of figure 4.3 when only the rigid 
body mode is included in the plant model (i.e. the flexible mode transfer functions 
are null), then 


Trpd ( - 5 ) 


since Cpd(s) 


a P r (s)C p d{s) a n T (s) 
~ 1 + PAsjC^s) d T (s) 

= p + ds 


q(p + ds) 

s 2 + s(b + ad) + ap 


(4.14) 

(4.15) 


where a>0, b > 0, p>0, d > 0 and a and b are taken from the plant 
analytical model. 

In order to configure T rp d(s) for a good step response we select the p and d 
gains as follows; let 

u 2 r = ap, (4.16) 

2 Cr^r = b+ad, (4.17) 

where u) T is the pole pair frequency and Cr their damping ratio. 

Fix £ r = 1 to ensure critical damping and select a set of increasingly large u T 
over a range of interest (e.g. 0.5 — + 4 Hz and converted to rad/ sec) and composed 
into a vector, 

U T 


= [ CcVj U r2 ... Um ] 


(4.18) 
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then equations 4.16 and 4.17 imply that 


•$3 

II 

(4.19) 

2( r U>ri b 

di = 

a 

(4.20) 

and p = | p] pj ••• Pn J * 

(4.21) 

and d = ^ d x d 2 ... d n j 

(4.22) 


where p and d are the n x 2 P and D monotonically increasing gains which 
correspond to u r . The higher the values of p, and d, the higher is u ri and the faster 
the closed loop step response. Note that the minimum u ri considered (i.e., for i = 1) 
is determined by equation ( 4.20) as 

^ (4 ' 23) 
in order to avoid any negative values for d{ since this may result an unstable 

Trpd(s). 


, 4.4.3 Root Locus Design 

The second step is to add the flexible modes to the plant model to produce 
P(s) as shown in figure 4.3 while maintaining the rigid body mode PD compensator 
class subject to the constraints of equations ( 4.19) and ( 4.20) so that the new 
torque to hub angle transfer function T p d{s ) is given by, 


w = 


P(s) 


nr{s) 


1 + Pi^C^s) d T (s ) 
and the PD compensator may be expressed as 


(4.24) 


u 

since yhub 
where Bj hift 


r - (p.Bl„ + d.B T )x 


(4.25) 

Bjhift 1 and = y/i’u6 = b t x 

(4.26) 

5(2) 5(3) ... 5(n) 0 


(4.27) 
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Figure 4.3: Hub Angle PD Feedback System 


In the following the first two flexible modes are considered and adding these 
extra plant states (four in all) will not destroy the closed loop system stability as 
the torque to hub velocity plaint transfer function remains strictly positive read. 

A root locus analysis of dr{s) yields the traice shown in figure 4.4 (the left plot 
is a more detailed view of the right plot close to the origin) where each open loop 
pole is marked with am x and the root loci emerge from these. In these root locii 
plots the proportioned and differential gains are varied in tandem over the elements 
of p and d but the underlying parameter being varied is the desired rigid body 
. mode closed loop pole pair frequency, ui r ,. The rigid body mode poles commence 
at s = 0 and s = —a and proceed to move adong the LHS real axis while the first 
flexible mode complex pole pair arc towards the real aixis until they meet, travel 
towards the origin and then separate at a lower frequency amd arc back towards the 
imaginary aixis. The second flexible mode complex pole padr move in a less draimatic 
but similar manner; the pole damping ratio increases to a maiximum (by a factor of 
three approximately) before decreasing as the poles return towards the imaginary 
axis with the pole frequency decreasing slowly. 
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Figure 4.4: Hub Angle PD Feedback Root Loci 


Table 4.1: PD Compensator Gains 


p 

d 

2.472 

1.287 


An indication of the flexible mode pole sensitivity to changes in the propor- 
tional and derivative gains is the density of root loci points which is very low in the 
case of the first flexible mode when the poles are on or near the real axis and for this 
reason the corresponding P and D gains are avoided. The final pole locations were 
selected as a compromise between maximizing the rigid body mode performance as 
indicated by the magnitude of the smallest real axis pole and minimizing the first 
flexible mode sensitivity to gain changes as well as maximizing its pole damping 
ratio (and hence minimizing its associated decay time constant). Both the open 
loop and closed loop pole parameters are given in table 4.2 (as well as those for the 
next controller which includes a second feedback loop) and the PD gains used are 
listed in table 4.1. 

The first flexible mode poles have a damping ratio close to 1/2, a natural 
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Table 4.2: Compensator Pole Parameters 


Mode # 

Jol 

fpd 

fapr 


Col 

Cpd 

Cspr 


Tot 

Tpd 

Tspr 

0 









0.92 

0.38 

0.40 

1 

2.83 

1.10 

1.06 


0.0229" 

0.81 

0.77 


2.5 

0.18 

0.19 

2 

7.15 

6.77 

6.71 


0.0115 

0.066 

0.133 


1.9 

0.35 

0.18 

3 

17.8 

17.8 

17.7 


0.0097 

0.0143 

0.029 


0.87 

0.63 

0.31 

4 

35.4 

35.4 

35.2 


0.0025 

0.0032 

0.0026 


1.7 

1.41 

1.74 

5 

58.0 

58.0 

57.9 


0.0016 

0.0017 

0.0013 


1.6 

1.62 

2.12 


1. subscript ol 
subscript pd 
subscript spr 

2 . t xx 


=> open loop plant parameter, 

=> PD compensated system parameter, 

=> PD + SPR compensated system parameter. 

_ i 
2 tt fxx^xx 


frequency of 1.1 Hz, and the slowest real axis pole has a time constant of 1.85 sec- 
onds.This closed loop complex pole pair will thus have some step response overshoot, 
but most of this will occur during the step rise time dictated by the slow real axis 
pole. 

The pole trajectories for the second flexible mode indicate relatively small gain 
sensitivity and less dramatic changes in damping ratios though for the selected p,- 
and d{ values their damping ratio has also increased by a factor of six approximately. 
Higher order modes show a progressively smaller increase in the damping ratios. 

The time domain state space representation, S pd for this closed loop system 
is as follows (assuming direct access to the states) 

x = ApdX + Bu (4-28) 

where A p d = A — B.p.Bj hi j t — B.d.B 7 (4-29) 

While A is in beam modal form (see equation ( 3.14)) A p d will have terms 
coupling the beam modes together and hence is not in beam modal form. 
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4.5 Synthesized Positive Real Output Controller (static feedback) 

The hub angle PD compensation feedback loop designed above is a form of full 
state feedback that provides good rigid body mode response with little or no over- 
shoot, approximately critically damps the first flexible mode and offers an infinite 
gain margin in the exact model case. However the higher order closed loop flexible 
modes remain relatively poorly damped (see table 4.1) and in order to add further 
damping to these flexible modes we can modify the PD full state feedback gains by 
closing a second feedback loop ; the following gain selection approach makes use of 
the time domain state space representation definition of strict positive realness in 
order to guarantee an infinite gain margin and a phase margin of [—90, 90] degrees 
for this second loop. In this example the exact gains are selected through trial and 
error, however there is a strong correlation between individual gain magnitudes (in 
a parametric form) and the corresponding closed loop flexible mode damping ratios 
which simplifies the selection process. 

• t A — 

Theorem la implies that open loop system P jpr (s) = B T PX p d(s ) (where 
Xpd(s) = (si — A p d)~ x B) will be a strictly positive real transfer function (and SISO 
since B is a column vector) if the square matrix P is the unique solution to 

A^P + PApd = ~L (4.30) 

where L > 0 and L* = L 

and the Passivity Theorem then guarantees that the close loop system T spr (s) = 
P apr (3)[I + P, pr (s)] -1 will be BIBO stable with an infinite gain margin. 

With regards to figure 4.5 X p d(s) is the torque reference (r(s)) to PD compen- 
sated plant state vector transfer function and Cspm = B T P is the output matrix 
synthesized specifically for this system so as to produce the SPR SISO transfer 
function P jpr (s), and since only static negative feedback is applied the correspond- 
ing closed loop system T spr (s) will have a phase margin of at least [—90, 90] degrees. 
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Figure 4.5: Hub Angle PD + SPR Feedback System 


The effect of this second loop is to apply negative feedback that can effect the damp- 
ing ratios associated with the new closed loop beam modes and we would now like 
to choose L so as to emphasize the poorly damped modes of through the output 
matrix C SVT \. However, since there exists coupling between beam modes (i.e. A p d 
is not in beam modal form due to the PD feedback terms) it is immediately clear 
how to do so. One solution is to convert T, p d to beam modal form using a similarity 
transformation, solve the Lyapunov matrix equation for a new positive real output 
matrix using this decoupled system and then transform the output matrix to suit 
the original state space. 

Let I = A p dX + Bu (4-31) 

and A T ?d P + PA pd = -L , l > 0 and I* = I, (4.32) 

where x = Mx , A pd = M~ l A^M , B = M~ l B and M is nonsingular : then 

B T P(sI-A pd )- 1 B = B T (M~ T PM~ l )(sI - Apd)~ l B (4.33) 

since y — B T Px = 1 3 t (M~ t PM~ 1 )x (4.34) 

and the RHS of equation ( 4.33) is a SPR transfer function since the Positive 
Realness Theorem guarantees that the LHS is a SPR transfer function. Hence we 
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can design a SPR output matrix B T P by applying the Lyapunov matrix equation 
to the PD compensated plant in a transformed state space representation and then 
derive the equivalent SPR output matrix C apr d = B T (M~ T PM~ l ) for the original 
state space representation. 

Now if M is chosen so that Apd is in beam modal form and L also restricted 
to the same block diagonal form then equation ( 4.31) shows that P will likewise be 
decoupled, i.e., if 



B 


CL\ 0 ... 0 

0 02 0 


X 


aili 0 ... 0 

0 02/2 ... 0 


0 0 ... a* 


0 0 ... a k h 


t>T h l 



a. , /, S S 2 * 2 , 6, € S 2 * 1 


(4.35) 


(4.36) 


then P = 


/, > 0, /* = and a, > 0 are scalars 


Pi 0 ... 0 

0 pj ... 0 


Pi € 3? 


2X2 


[0 0 ... pk J 

where af + a^ai = — a;/, 

and B T {M- T PM- x ){sI-A pd y 1 B 


(4.37) 


(4.38) 


= ^ ctibj p,(s/ - a,) 1 6, (4.39) 

1=1 

and i € {0, ..., &i}. The RHS of equation ( 4.39) shows that the the contribu- 
tion to the s.p.r. transfer function B T (M~ T P — A p d)~ x B by the i tfl beam 
mode SPR transfer function 6fp,(s/ - a,) -1 ^ is scaled by the positive constant 
a,- which may be selected as desired while the form of /,• determines the relative 
weighting each of the two states associated with the same beam mode. 
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Table 4.3: Lyapunov Equation relative Weights 



oro/ao 

Ct\/oc 0 

aj/ao 

Q^/oro 

Coupled L 

1 

20 

60 

80 

Decoupled L 

1 

1 

26 

6 


We now have two new compensators 


u 

= 

r - ( p.Bj hi j t + d.B T + C sp r)x 

(4.40) 

cv 

A 

B t P 

(4.41) 

and u 

=r 

r - ( p.Bj hi j t + d.B T + CsPRd)x 

(4.42) 

Capri 

A 

B t (M- t PM~ 1 ) 

(4.43) 


where P is the solution to the Lyapunov matrix 4.30 and P is the solution 
to the decoupled form ( 4.31) and an example of each is selected in the following 
paragraphs. 

4.5.1 Selecting Csprd (decoupled version) 

In order to emphasize the third and fourth PD closed loop beam modes several 
L matrices were tested with U the 3f£ 2 identity matrix and 03 and oc± large w.r.t. cx 1 
and c*2. Increasing for example 03 increased the feedback gains associated with the 
third beam mode and the corresponding closed loop mode damping ratio ; in this 
way the set of a, shown in table 4.3 were selected so as to increase the closed loop 
damping ratios of the third and fourth beam modes by a factor of approximately 
two over those of the PD compensated case while only slightly altering those of the 
first and second beam modes. The new closed loop pole parameters are not shown 
but are very similar to those given in table 4.2 for the PD+SPR case. 
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4.5.2 Selecting C apr (coupled version) 

For the SPR output design based upon the original coupled system Spo the 
form of L chosen was that of L in the previous example (i.e., as in equation ( 4.35) 
with U the 3J 2 identity matrix) with again small values for qj and a -2 and large 
values for a 3 and a 4 . Although there is no longer an exact correlation between a l 
and the contribution of each beam mode of S pd to the synthesized SPR output 
there is an approximate effect if the intermodal coupling is small as is the case for 
the third and higher order flexible beam modes. Again a 3 and a 4 were selected so 
as to increase the closed loop damping ratios of the third and fourth beam modes 
by a factor of approximately two over those of the PD compensated case while only 
slightly altering those of the first and second beam modes and the results are shown 
in table 4.2 while the relative values for a, are given in table 4.3 

4.6 State Estimation 

The experimental testbed currently contains a hub angle sensor and four beam 
strain gages : any other output such as a velocity signal or a synthesized passive 
output must be obtained by direct filtering of the existing outputs or reconstructed 
from an observer estimate of the plant states. The latter was chosen so as to 
provide a complete plant state estimate from which any desired output may be 
generated. Under ideal conditions a close match may be obtained between the 
actual and estimated state within the bandwidth of the observer, however modeling 
errors and unmodeled features can distort the estimated state and in this section 
we look at such effects due to severed uncompensated beam modes (i.e., some higher 
order beam modes not modeled at the compensator design stage) in the presence of 
three different types of closed loop observers. In particular, we apply PD negative 
feedback based upon the estimated state of k\ beam modes and then examine the 
transfer function from the single torque reference signal to the single synthesized 
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SPR output (as designed in the previous section) in the presence of k 2 extra beam 
modes not modeled when the observer and compensation were designed. The three 
observer designs used are 1. full state pole placement, 2. full state LQE and 3. 
reduced order Luenberger. 

An open loop observer using an inverted plant output matrix based upon 
existing sensors cannot reconstruct plant velocity states since they are decoupled 
from these output while bandpass or highpass filtering of the existing outputs require 
relatively fast poles (in order to minimize the associated phase shifts) and hence will 
amplify high frequency noise. A closed loop observer however may reconstruct the 
plant velocity states (as well as amplitude states) with poles placed close to the 
highest modal frequencies of interest while maintaining a low phase estimate error. 

In the following the observers and compensators will be designed using a beam 
model of k x beam modes and then tested with a beam model of k x + k 2 beam modes 
(see figure 4.6) and therefore it will be useful to partition the plant model as follows; 


Xi 



A x 0 


Xx 

+ 

Bx 

*2 _ 


° a 2 


. x 2 . 


b 2 

y = 

9 


Xx 





. j 


x 2 




(4.44) 


(4.45) 


where A x € % 2k ' x2k ', A 2 € % 2k * x2k *, B x € % 2k ' xl , B 2 e 3 * 2fcjXl , C x € & px2 *-, 
C 2 € 3J px2fc *, u € & is the torque input signal and y € 3£ p is the plant output vector: 
[Ai, B\, C\ ] is the k x beam mode compensated system while [A 2 , B 2 , C 2 \ represents 
the extra beam modes unmodeled at the observer/compensator design stage. 

The PD negative feedback is a form of full state feedback, i.e. 

u = r — Kpdii (4-46) 

where K p d = p.Bj hift + d.B T ( see equation ( 4.27) ), (4.47) 


x x is the observer estimate of x x and r is the torque reference signal. The 
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Figure 4.6: SPR Open Loop System with Unmodeled Modes 


synthesized positive real output is given by 

tispr = r xi ■ •• (4.48) 

where y, pr € 3? and C S pr € 3ft lx2 *‘. 

4.6.1 Full State Observers with Uncompensated Modes Effects 
The full state closed loop observers tested may be represented by 

ii = Aiii + Biu + L(y - y) (4.49) 

y = Cifi (4.50) 

where L € 3ft 2 * 1 xp : L will be selected using pole placement (e.g., L pp ) and also 
LQE (e.g., Li qe ) techniques. 

We can examine the combined effects of the unmodeled plant beam modes and 
the observer on the transfer function by studying the estimated state error 

dynamics where we define the 3ft 2 * 1 xl state error vector as 


e 


X\ — X\ 


(4.51) 




55 


Equations ( 4.51), ( 4.44) and ( 4.46) imply that 

y pr (s ) = Ospri i( 5 ) = C 3 prTi p d(s')r(s) + C $ pr 'I , e i(s')t(3') (4.52) 

where Tipd(s) = (si — A\ p d) l B\ , (4.53) 

T ti (s) = I + {sl — Ai p d)~ l BiKpd, (4.54) 

and A\pd = A\ — B\Kpd ■ (4.55) 

Ti p d(s) is the designed PD closed loop transfer function (torque reference to 
modeled plant state) while T em (s)e{s) represents the associated (e)rror dynamics for 
the (i)deal case which in the absence of unmodeled beam modes will be exponentially 
stable and decoupled from v since equations ( 4.49) and ( 4.51) show that 

e = A 0 e (4.56) 

where A 0 = A\ + LC\. (4-57) 

Now when unmodeled beam modes are included in the plant e is driven by the 
extra plant states X2 and the corresponding error equation is given by 

e = A 0 e-LC 2 x 2 (4.58) 

or e(s) = — T ue /(s)x 2 (s), (4.59) 

where T ue /{s) = (si — A 0 )~ l LC 2 . (4.60) 

(4.61) 

Tue f is the transfer function between the (u)nmodeled plant states and the 
modelled state (e)rror vector for the (f)ull state observer case. 

Lastly, combining equations ( 4.59), ( 4.46) and ( 4.52) yields 

CsprXiis) = C spr T lp d(s)(I + (I- UsWpdr'T.WKpdHs) (4.62) 

where T a (s ) = T e ,(3)r ue /(s)T2p < i(s) (4.63) 

and T 2p d(s) is the unmodeled state counterpart to Ti p d(s) 
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i.e. T 2p d{s ) — {si ~ A 2p d) 1 B 2 (4.64) 

where A 2p d = A 2 — B 2 K P d (4.65) 

Equation ( 4.62) shows the effect of the unmodeled dynamics (i.e., T 2p d(s)) on 
the synthesized positive real output C^x^s) (which of course may no longer be 
strictly positive real): if T 2p d(s) is null then this equation reduces to the ideal case 
(i.e equation 4.52 with zero error initial conditions). It is desirable to maintain the 
positive realness of C 3pr xi{s ) in order to preserve the associated obtainable gain and 
phase margins when this second feedback loop is closed and hence equation ( 4.62) 
suggests that L and C 3pr be selected so as to maximize 


inf Re(C, pr T lp d(ju)[I + {I- T a {ju)K pd )~'T a {juj)K pd \r{ju)) (4.66) 


(as well as to meet other observer design criteria and emphasize selected beam 
modes). 

4.6.2 Reduced Order (Luenberger) Observer with Uncompensated Modes 
Effects 

The reduced order (Luenberger) observer may be represented by 

zi = A ro zi + Eu + Gy (4-67) 

(4.68) 

(4.69) 

(4.70) 



_ _ 

-1 

r 1 


c, 


y 

x x = 





r 


z 


where A ro T = TAi — GCi, 
E = TBi 


and T is also constrained such that 


C x 

T 


is nonsingular. 


The analysis of the previous section may be carried out in order to determine 
the transfer function relating the torque reference r(s) to the synthesized SPR output 
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C7 apr x 1 (s) with only one change, namely the relationship between x 2 and e. Unlike 
the full state observer case the term xi(s) does not appear explicitly in the reduced 
order observer state equations so we work with the linear combinations that do, 
namely CiZi(s) ( = y ) and Txi(s) ( = z ) and relate these to e and then equation 
( 4.68) gives the desired transfer function. 


Ciii = y = C\Xi + <7 2 x 2 (4.71) 

=* Cje = C 2 x 2 (4.72) 

and equations ( 4.67), ( 4.69), ( 4.70) and 4.51 imply that 

Te = z-TZx = A ro Te + GC 2 x 2 (4.73) . 

Te(s) = (si - Aro)- l GC 2 x 2 (4.74) 


Combining ( 4.51), ( 4.68), ( 4.72) and ( 4.74) gives 

I a r 1 ' 

e(s) = 


Cl 

T 


Cie(s) 

Te(s) 


(4.75) 


-1 


Cx 

T 


(sI-A ro )~ l G 


C 2 x 2 (s) = T uer {s)x 2 (s). (4.76) 


T uer (s) (r for (r)educed order Luenberger observer) may now replace T ue j(s) 
in equation ( 4.62) above to give the new r(s) to C^Xi(s) transfer fimction when 
the unmodeled beam modes axe included. 


CHAPTER 5 

Simulation and Experimental Results 

In the previous chapter two controllers were developed for the analytical linear model 
of the flexible beam, PD and PD plus synthesized strictly positive real feedback, and 
in the following section the step responses of the corresponding closed loop systems 
are examined using both linear model simulations and experimental tests on the 
actual testbed (and therefore in the presence of unmodeled linear and non-linear 
dynamics and noise). Also, distortion effects under the three closed loop observers 
are evaluated and the best observer selected. 

The compensators and observers are converted to discrete form using the bilin- 
ear Tustin transformation and a sampling frequency of 500 Hz which is over twenty 
five times higher than the highest compensated beam modal frequency : the same 
descretized models are used for both computer simulations and experimental tests. 
In each case reference trajectories are applied in the form of open loop plant modal 
states. 

In each case the closed loop system structure is very similar to that of the open 
loop plant; e.g., the hub angle PD and passivity based feedback systems contain a 
rigid body mode (though there is no longer a pole at the origin) and flexible modes 
with new center frequencies and damping ratios. 

5.1 Simulation Results with Ideal Observer 

The PD full state controller designed in the previous chapter provides a good 
rigid body mode response and large damping for the Second closed loop flexible 
mode. However higher order closed loop flexible beam modes remain poorly damped 
and a second feedback loop based upon a synthesized strictly positive real output is 
used to increase the corresponding damping ratios by a factor of two approximately 
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Figure 5.1: Experimental Controller Architecture 

for the second and third flexible beam modes. In the absence of observer effects the 
differences between the two examples of these SPR based feedback are negligible 
and hence only the first version will be considered in this section. 

Figure 5.2 shows the hub angle (a), first strain gage (b) and tip displacement 
(c) step responses (simulated) for the two controllers (PD and PD + S.P.R.) with a 
reference step size of 45 degrees. The hub angle and tip displacement 90% rise times 
are similar for each controller (approximately 1.8 seconds) since these are dominated 
by the closed loop rigid body and first flexible beam modes which are similar for 
both controllers. An increase in damping for the second and third flexible beam 
modes is evident in the strain gage plot for the PD + S.P.R controlled system which 
shows less residual oscillations. Figure 5.3 is another view of the first strain gage 
(a) and tip responses (b) starting at t = 4 seconds that more clearly shows residual 
oscillations due to closed loop flexible modes : again the PD + S.P.R controller 
results in greater damping. 




angle <deg> 
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5.2 Closed Loop Observers Effects 

When observer and unmodeled beam modal dynamics are added to the com- 
puter simulations the actual transfer function from torque reference to synthesized 
SPR output, fipd(s) = C spr ii{s)/r(s) (see equation ( 4.62)), differs from the ideal 
one, Tipd(s) (see equation ( 4.53)) and these non- ideal transfer functions are com- 
pared in figures 5.4, 5.5 and 5.6 to Ti p d(s). For each of the three closed loop 
observers selected the error ti p d{s) - T lp d{s ) is large in the bandwidth of the un- 
modeled modes and a large phase shift (much greater than —90 degrees) also occurs 
which means that Ti p d(s) is no longer SPR: in fact phase shifts greater than -180 
degrees occur with gain magnitudes larger than unity and hence the stability of the 
closed loop system is no longer guaranteed. 

Within the bandwidth of the compensated plant each observer results in a 
different transfer function error with the reduced order Luenberger model providing 
the best match, the LQE designed full state observer the worst match and the pole 
placement full state observer coming somewhere in between. 

Figure 5.7 shows the experimental hub angle and first strain gage step re- 
sponses for each of the three observer cases tested (using a PD + SPR controller) 
and the reduced order Luenberger and full state pole placement models produce 
the best results (they are almost identical to each other). The hub reference plot 
is shown as a dashed line. Because of the good results obtained with the reduced 
order Luenberger observer both in simulation and experiment it is the only one used 
in the next section when comparing controllers, however it is emphasized that each 
observer was designed through trial and error and hence we cannot conclude from 
these results that one observer type is categorically ’better’ than another in the 
presence of unmodeled beam modes and other model mismatches. 

The transfer function error T\pd(s) — T lp d(s) also depends on the selected syn- 
thesized SPR output and figure 5.8 compares Ti p d(s) to TU*) for two different 


Frcq <H z>. 


Figure 5.5: SPR T.F. Bode Plots with Pole Placement F.S.Observer 
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Freq <Hz>. 


Figure 5.6: SPR T.F. Bode Plots with LQE F.S. Observer 

such outputs; one design uses the PD closed loop plant in coupled form while the 
other used the decoupled system and both simulations are with the reduced order 
Luenberger observer. In each case the transfer function gain error within the band- 
width of the compensated beam modes is small however the decoupled model based 
SPR output design contains a non-minimum phase zero that increases the phase 
shift by —180 degrees approximately with respect to the minimum phase case. For 
this reason the coupled model based SPR output was used in all other simulations 
and experimental tests although it is again emphasized that trial and error was used 
in the design process, in this case when selecting the synthesized SPR output. 

5.3 Experimental Results with Reduced Order Luenberger Observer 

Figures 5.9 and 5.10 show the hub angle (a) and first strain gage (b) step re- 
sponses (experimental) for the two controllers (PD and PD + SPR) with a reference 
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Figure 5.10: Step Response Experimental Results Closeup 

step size of 45 degrees. The results are similar to those of the computer simulations 
which is encouraging however residual oscillations are larger and a steady state error 
exists due to friction effects at the hub. 

As in simulations the PD + S.P.R. controlled system has less residual oscil- 
lation than the PD case and the an estimate of the increase in damping ratios of 
closed loop poles for the second and third flexible modes may be obtained by ex- 
amining magnitude Bode plots derived from experimental results when the beam 
was excited simultaneously at a range of frequencies, see figure 5.11. In the case 
of an isolated flexible mode the peak transfer function Bode gain under closed loop 
differential feedback is inversely proportional to the closed loop pole damping ratio. 
This effect is distorted with the introduction of other transfer function dynamics, 
e.g., other flexible modes, however if these modes are sufficiently well spaced apart 
in the frequency domain and lighlty damped than the transfer function close to 




67 


3 


.s 



frequency <Hz> 



Figure 5.11: Experimental Bode Plots 


the center frequency of an individual flexible mode closely approximates that of the 
corresponding mode and hence in figure 5.11 (a) the reduction in gain magnitudes 
between the two plots at the flexible mode center frequencies are an approximate 
indication of the increase in close loop pole damping ratios: an attenuation by a 
factor of approximately two for the second and third closed loop flexible modes is 
as predicted by the simulations. Figure 5.11 (b) is a log scale version of figure 5.11 
(a) that clearly shows the zeros which are common to the open loop system (see 
e.g., figure 3.7 part (c)). 




CONCLUSIONS 


The hub angle PD compensator derived from the linear normal mode model produces 
very good results in terms of the hub angle step response which has no overshoot 
and little residual oscillation in both simulations and real tests The same is also 
true for the beam tip displacement for the simulated case and although this cannot 
yet be directly verified for the real testbed case low residual oscillations in the tip 
displacement are implied by the low residual oscillations in the strain signals. 

Adding a second SPR static gain full state feedback loop can improve the step 
response performance further in a robust manner by increasing the flexible mode 
pole damping ratios. Such compensation has a phase margin of at least [—90, 90] 
degrees in the ideal case and this is found to be also true for simulations that 
include observer effects and some unmodeled modes within the bandwidth of the 
compensated modes though not at higher frequencies. This suggests the addition 
of a low pass filter in the SPR feedback loop that maintains a low phase shift and 
nearly constant gain in the compensated mode bandwidth in order to maintain the 
low frequency closed loop performance already achieved, with a gain roll off at higher 
frequencies to attenuate unwanted effects. 

Of the three closed loop observers tested the two with poles placed above the 
highest compensated mode center frequency produced the least distortion on the 
synthesized SPR transfer function and the best experimental closed loop perfor- 
mance. 

The controller design procedure used appears basically sound and produced 
good results however it is fractured and involved trial and error selections. The task 
now is to synthesize these observer /compensator blocks either separately or together 
in a systematic and optimum fashion. 
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APPENDIX A 

Beam and Actuator Physical Parameters 


This appendix contains four tables which list the physical parameters associated 
with the flexible beam testbed and most are used when deriving the beam model in 
chapter 3. 


Table A.l: Aluminum Beam Material Parameters 


Parameter 

Symbol 

Units 

Value 

Width 

w 

m 

1.5875e-3 

Height 

h 

m 

1.03e-l 

Length 

l 

m 

1.098 

Density 

P 

Kg/m 3 

2.713e3 

Youngs Modulus 

E 

Kg/m 2 

6.895el0 


Table A. 2: Beam Derived Parameters 


Parameter 

Symbol 

Units 

Value 

Formula 

Inertia (about hub axis) 

h 

Kg.m 2 

1.197e3 

p ‘ 73 

Cross Sectional Area 

A^ C a 

m 2 

1.6345e-4 

w.h 

Mass per Unit length 

P r 

Kg/m 

4.4327e-l 

p'Abca 

C.S. Area Inertia 

la 

m 4 

3.4316e-ll 

A hca .w 2 j 12 

Mass 

M b 

Kg 

4.436e-l 

w.h.l.p 


Note that the actual hub inertia value used in the model is 2.5 x 10 -2 Kg.m 2 
(rather than the conservative 2.363 x 10 -2 Kg.m 2 value computed in table A.4) as 
this was the best approximation at the time. Each mass and inertia was calculated 
from the dimensions and properties of the system components is only an approximate 
value that does not take into account fine detail or small objects such as screws etc. 
Also, the extra mass and inertia due to the light sensor detector and optics currently 
mounted at the hub are not included in table A.4. 
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Table A.3: Motor Parameters 


Parameter 

Symbol 

Units 

Value 

Shaft Inertia 

Im 

Kg.m 2 

1.201e-4 

Shaft Mass 

M m 

Kg 

— 

Torque Constant 

KT 

Nm/A 

2.370e-l 

Continuous Stall Torque 

TS 

Nm 

1.2407 

Peak Torque 

TP 

Nm 

1.358el 

Electrical Time Constant 

TE 

sec 

1.4e-4 


Table A. 4: Hub Derived Parameters 


Item 

Symbol 

Mass 

M h 

Kg 

Inertia 

h 

Kg.m 2 

# 

Motor Shaft 


1.201e-4 

1 

Hub Collar 

4.798e-2 

3.869 e-6 

1 

Encoder Collar 

— 

— 

1 

Bellows 

— 

— 

1 

Potentiometer 

— 

— 

1 

Hub Base 

1.126 

6.707e-3 

1 

Hub Base Strengthener 

— 

— 

1 

Counter Weight 

7.476e-l 

1.482e-2 

1 

Counter Weight Bar 

4.123e-2 

7.047e-4 

1 

Balance weight Bar 

2.61 le-2 

3.263e-4 

2 

Weight Bar Block 

3.619e-3 

3.000e-5 

3 

Beam Post 

3.941e-l 

9.214e-4 

2 

Hub Assembly total 

2.386 

2.363e-2 
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